Optical properties of self-organized wurtzite InN/GaN quantum dots: 
A combined atomistic tight-binding and full configuration interaction calculation 
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In this work we investigate the electronic and optical properties of self-assembled InN/GaN quan- 
tum dots. The one-particle states of the low-dimensional heterostructures are provided by a tight- 
binding model that fully includes the wurtzite crystal structure on an atomistic level. Optical dipole 
and Coulomb matrix elements are calculated from these one-particle wave functions and serve as an 
input for full configuration interaction calculations. We present multi-exciton emission spectra and 
discuss in detail how Coulomb correlations and oscillator strengths are changed by the piezoelectric 
fields present in the structure. Vanishing exciton and biexciton ground state emission for small 
lens-shaped dots is predicted. 

PACS numbers: 78.67.Hc, 73.22.Dj, 71.35.-y 



In recent years, semiconductor quantum dots (QDs) 
have been the subject of intense experimental and the- 
oretical research. As a new material system, group-Ill 
nitride based devices are of particular interest due to 
their wide range of emission frequencies from red to ul- 
traviolet and their potential for high-power electronic 
applications 1 - 2 . Being a technologically promising sys- 
tem, we study self-assembled InN/GaN QDs, which are 
typically grown by molecular beam epitaxy in Stranski- 
Krastanov growth mode. A theoretical description of 
the one-particle states in terms of a tight-binding (TB) 
model is presented, which provides a powerful approach 
to the electronic states of low-dimensional heterostruc- 
tures on an atomistic level 3 ^. For the calculation of op- 
tical absorption and emission spectra, full configuration- 
interaction (FCI) calculations^ are used to obtain a 
consistent description of correlated many-particle states. 
The calculation of dipole and Coulomb matrix elements 
from the TB one-particle wave functions facilitates the 
combination of these two approaches and allows us to 
investigate optical transitions between the interacting 
many-particle states of a QD with parameters obtained 
from a microscopic model. For the investigated small 
lens-shaped InN/GaN QDs, we report a negligible ex- 
citon and biexciton ground state emission whereas at 
higher excitation conditions strong emission from three 
to six exciton complexes is obtained. 

We consider lens-shaped InN QDs, grown in (000 In- 
direction on top of an InN wetting layer (WL) and 
embedded in a GaN matrix. Their circular symme- 
try around the z-axis (diameter d = 4.5 nm, height 
h = 1.6 nm) preserves the intrinsic C^ v symmetry of 
the wurtzite crystal. For the WL we assume a thick- 
ness of one lattice constant. We apply a TB-model 
with an sp 3 basis |a,R), i.e., one s-state (a = s) and 
three p-states (a = PxiPyiPz) per spin direction at each 
atom site R. In contrast to most other III-V and II- VI 
semiconductors, one can neglect spin-orbit coupling and 
crystal-field splitting in InN and GaN2*£. We include 



non-diagonal elements of the TB-Hamiltonian matrix up 
to nearest neighbors and use the two-center approxima- 
tion of Slater and Koster— which yields 9 independent 
TB-parameters. These parameters are empirically deter- 
mined such that the characteristic properties of the bulk 
bandstructur o 9 : 10 in the vicinity of the T point are repro- 
duced. With these TB-parameters, the QD is modeled 
on an atomistic level where the parameters for each site 
are set according to the occupying atoms (N, In, Ga). 
At the InN/GaN interfaces averages of the parameters 
arc used and effects of the surfaces in the finite-size su- 
percell are removecUi. The spontaneous polarization in 
the wurtzite crystal structure lies within growth direc- 
tion: P = Pe z . Additionally, a strain-induced piezoelec- 
tric field occurs that is quite strong for the investigated 
InN/GaN heterostructures. The piezoelectric field is de- 
termined by solution of the Poisson equation. The strain 
contribution to the polarization is approximated in a way 
following Ref. 0. For our chosen dot geometry, even a 
more sophisticated inclusion of strain effects^ will gen- 
erate merely small lateral contributions to the piezoelec- 
tric fieldi^, which are therefore neglected in the follow- 
ing. Lattice mismatch parameters and strain tensors are 
taken from Ref.ll5l the small thermal strain contribution 
is neglected^. The calculation yields a reasonable value 
of 5.5 MV/cm for the electric field inside the QD. The re- 
sulting electrostatic potential is included in the TB model 
as a site-diagonal potential V p {v) = —e (j)p{r). This 
method has successfully been applied to quantum welli& 
and QDi4 structures before. By including the piezoelec- 
tric field the quantum confined Stark effect (QCSE) and 
its influence on the Coulomb matrix elements and the 
oscillator strengths can be studied. 

The discussed QD confines three bound states for the 
electrons. They are included in the FCI calculations, to- 
gether with the three lowest one-particle hole states that 
are spectrally well separated from the other localized hole 
states. These three one-particle states for the electrons 
and holes and their energies are depicted in Fig. ^ As 
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FIG. 1: (Color online). The QD geometry is shown from 
atop. The structure is visualized and isosurfaces of the charge 
density for the three energetically lowest electron (left) and 
hole (right) states are included for 10 % (blue) and 50 % (red) 
of the maximum value. For the holes the atomistic structure 
of the wurtzite crystal becomes most apparent for the excited 
p-states. The corresponding energies (E^' 2 3 ) of electron and 
hole states measured from the valence band maximum of bulk 
GaN and the atomic orbital character for each wave function 
are given. The dominant contributions are highlighted. 



the lowest state for electrons and holes is invariant un- 
der rotation by 2ir/3, it is denoted as s-state. The two 
excited states are classified as p-states according to their 
symmetry properties and their two-fold degeneracy. The 
hole states show a strong band mixing visible in the or- 
bital character (Fig.^), which is in agreement with other 
multi-band approaches'^. 

As emphasized above, a TB model represents an atom- 
istic approach to describe the electronic structure of low- 
dimensional heterostructures. However, explicit knowl- 
edge about a basis set of localized Wannier states is 
not required for the calculation of one-particle energies 
and wave functions. Only the basic assumptions about 
the localized (atomic) orbitals, i.e. symmetry, spatial 
orientation^, and orthogonality, enter the TB Hamilto- 
nian. Nevertheless, for the calculation of dipole and 
Coulomb matrix elements one needs - in principle - the 
localized basis states. For the Coulomb matrix elements, 
however, which are dominated by long-range contribu- 
tions, the explicit knowledge of the atomic orbitals is in 
practice not required. This is because the structure of 
the localized orbitals is of significance only for on-site 
and nearest-neighbor interactions, which contribute less 
than 5% to the total Coulomb matrix elements. These 
findings are in agreement with Ref. Thus, the ma- 
trix elements are approximated by a sum over the TB 
coefficients at atom sites R, R' with orbital indices a, (3: 
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The labels i,j,k,l refer either to electron or to hole states 
in case of the repulsive electron-electron and hole-hole 
interaction, or i, I label electron and j, k hole states for 
the attractive electron-hole interaction. The considerably 
smaller matrix elements of the electron- hole exchange in- 
teraction are neglected. The electronic charge and the 
vacuum dielectric constant are denoted by eo and eo, re- 
spectively. We use the InN dielectric constant s r = 8.4 
according to Ref. since the wave functions are almost 
completely confined inside the QD. 

For the calculation of dipole matrix elements d^j 1 = 
eo(ipf \r\il>j), the explicit structure of the localized or- 
bitals is required as the dipole-operator has mainly lo- 
cal character. Standard Slater orbitals^ have been used 
in earlier calculations^ within orthogonal TB models. 
While they include the correct symmetry properties, the 
missing orthogonality limits their applicability. To over- 
come this problem, we use numerically orthogonalizcd 
Slater orbitals. To properly treat the slight non- locality 
of the dipole operator^ and the anion-cation structure of 
the crystal, the matrix elements are calculated including 
up to second nearest neighbors. The only relevant dipole 
matrix elements are edli 1 = edf„ and edfr . = edr*,,, 

SJJx ° Py P^e * Py » ' 

where e = 1/V2(1, 1,0) denotes the light polarization 
vector. All other matrix elements vanish due to the over- 
all symmetry of the connected one-particle statesSSiS. 
The resulting optical selection rules are in strong con- 
trast to what is known from many other III-V and II- VI 
heterostructures and cannot be explained within an one- 
component effective-mass approach&Si. 

The single-particle states and Coulomb interaction ma- 
trix elements are used in the FCI calculations to deter- 
mine the multi-exciton states. In a second step, Fermi's 
golden rule is evaluated for dipole transitions between 
these Coulomb-correlated states in order to obtain the 
multi-exciton emission spectra^. The results for an ini- 
tial filling of the dot with one up to six excitons are de- 
picted in Fig. with (solid line) and without (dotted 
line) the piezoelectric field. In the emission spectra two 
clusters of peaks are clearly visible, one on the high en- 
ergy side, fojj > 1.24 eV, and one on the low energy side, 
hw < 1.3 eV (explicit numbers given in the text refer to 
the results including the piezoelectric field). As a char- 
acteristic feature, the high (low) energy transitions orig- 
inate from a recombination process involving an electron 
(hole) in an excited state and a hole (electron) in the 
ground state. As the dipole matrix element d e Jl is neg- 
ligible and the exciton and biexciton ground states are 
dominated by configurations where all the carriers are in 
the s-shell, the corresponding transitions IX — > OA and 
2A — > 1A remain dark. The gap between the two sets of 
clusters is, to a first approximation, given by the differ- 
ence in the involved one-particle energies, which is then 
renormalized by the Coulomb interaction. In a free parti- 
cle picture each of the two clusters would collapse to one 
line. The splitting within the clusters can be attributed 
to transitions with different spin configurations of the 
final states. These configurations are energetically sepa- 
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FIG. 2: (Color online). Emission spectra for a quantum dot 
with different number of excitons, with (solid line) and with- 
out (dotted line) piezoelectric field. Initial states are ground 
states with total spin z-component S z = 0. 



rated by approximate integer multiples of the Coulomb 
exchange integrals V?Jr ee a „ ~ 18 meV for the low en- 

&Psc , y **Px ,y 

ergy cluster, and Vj^ hh spn; ss 4meV for the high energy 
cluster. For clarification let us consider the 3X — > 2X 
transition. The 3X ground state is dominated by config- 
urations with both electrons and holes having two car- 
riers in the s- and one in the p-shell. Recombination of 
the p-electron (hole) with the s-hole (electron) leaves be- 
hind two electrons (holes) in the s-shell as well as one 
hole (electron) in the s-shell and one in the p-shell. The 



latter can either form a singlet or a triplet state, which 
leads to a splitting of around 8meV (36meV), and is in 
accordance with the FCI-calculation. As the number of 
e-h-pairs is increased one observes an overall blue shift 
of the transitions. This shift is less pronounced with- 
out the piezoelectric field and can be explained in terms 
of the Hartree-Fock contributions to the interaction. The 
in general somewhat larger peak height of the low-energy 
transition is explained by the fact that the involved dipole 
matrix element d e Jl is larger than df! 1 „, which enters 

&Px ,y Px ,y& 

the high energy transitions. 

The inclusion of the piezoelectric field for the strained 
wurtzite crystal structure gives rise to a QCSE which cre- 
ates a strong (about 220 meV) redshift of the one-particle 
gap energy. Additionally, the Coulomb matrix elements 
are strongly modified and the oscillator strengths are ap- 
proximately halfed due to the spatial separation of elec- 
tron and hole wave functions. 

In conclusion, we successfully combined two state-of- 
the-art approaches, the atomistic tight-binding model 
and full configuration-interaction calculations, to inves- 
tigate the optical properties of the technologically very 
promising InN/GaN QD system. Multi-exciton emission 
spectra are calculated with microscopically determined 
input parameters, which reveal the strong influence of 
bandmixing effects on the optical transitions between the 
Coulomb correlated many-particle states. As an impor- 
tant consequence for future optoelectronic applications 
we predict vanishing exciton and biexciton ground state 
emission for small lens-shaped InN/GaN QDs. 
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